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ABSTRACT 

We report new 2 cm VLBA images of the inner radio jet of M87 showing a limb brightened structure and 
unambiguous evidence for a faint 3 mas long counter-feature which also appears limb brightened. Multi-epoch 
observations of seven separate jet features show typical speeds of less than a few percent of the speed of light, 
despite the highly asymmetric jet structure and the implications of the canonical relativistic beaming scenario. 
The observed morphology is consistent with a two stream spine-sheath velocity gradient across the jet, as might 
be expected from the recently discovered strong and variable Te V emission as well as from numerical modeling 
of relativistic jets. Considering the large jet to counter-jet flux density ratio and lack of observed fast motion in 
the jet, we conclude that either the inner part of the M87 jet is intrinsically asymmetric or that the bulk plasma 
flow speed is much greater than any propagation of shocks or other pattern motions. 

Subject headings: galaxies: active — galaxies: jets — galaxies: individual (M87) — radio continuum: galaxies 
— acceleration of particles 



1. INTRODUCTION 

The peculiar galaxy M87 in the Virgo cluster was among 
the first to be recognized as a powerful source of radio emis- 
sion. More than 40 years ago, Sh klovskvl (Q~964) argued that 
since radio galaxies such as Cygnus A typically have symmet- 
ric lobes, the jets which feed the radio lobes are likely to be 
intrinsically two-sided, but that they would appear one-sided 
due to differential relativistic Doppler beaming. This simple 
picture is now widely accepted as the basis of unified models, 
which try to explain many of the different properties of AGN 
as due to relativistic beaming a nd projection effects, ra ther 
than to intrinsic differences (e.g.. lUrry & Padovanil 19951) . 

M87 remains of great interest today, since there is strong 
observational evidence for a 3 x 10 9 M Q black hole located 
at the galactic nuc l eus thought to power the relativistic jet 
(lHarms et al.lll994t iMacchetto et all 1 19971) . Moreover, at a 
distance of only 16 Mpc, (1 mas = 0.08 pc; 1 mas per year 
= 0.25c) M87 is one of the closest radio galaxies, and as such 
it has one of the few jets which can be well-resolved on sub- 
parsec scales in a direction transverse to the flow. 

The M87 jet has been studied over a wide range of 
the electromagnetic spectrum, incl uding imaging w ith HST 
dBiretta et al.ll 19991) and the VLA jOwen et alii 19891) as wel l 
as with Chandra dWilson & YangN2002t lHarris et al.| [2006). 
The observed morphology at radio, optical, and even X-ray 
wavelengths appears very similar, suggesting a common syn- 
chrotron radiation mechanism at all wavelengt hs and a com- 
mon s pectral index of -0.67 throughout the jet (Perl man et al.l 
1200 ll) . At least at optical and X-ray wavelengths, the elec- 
tron lifetime down the jet is much shorter than the travel time 
from the nucleus, suggesting conti nual in situ accelera tion of 
relativistic particles within the jet (Perl man et al.l l2001). 
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In this paper we report on observations made with the 
NRAO Very Long Baseline Array (VLBA) at 2 cm wave- 
length. High dynamic range images constructed from obser- 
vations made in the year 2000 describe the two-dimensional 
structure of the jet out to nearly 0.2 arcsec (16 pc) and show 
the presence of a faint counter-feature. These observations 
were complemented by regular observations made with lower 
sensitivity between 1995 and 2007 to study the outward flow 
within the inner part of the radio jet. 

2. THE OBSERVATIONS AND IMAGING 

M87 has been regularly observed with the VLBA 
since 1995 as part of the 2 cm VLBA survey 
dKellermann et alj 119981) and the more recent MOJAVE 
program dLister & H oman 2005). In these programs, at each 
epoch we observed each source for a total of about one 
hour, with multiple observations spaced over a wide range 
of hour angle. For M87, we have obtained a total of 21 
images between 1995 and 2007. Typically the rms noise in 
each image is about 0.3 mJybeam" 1 . A s d escribed in more 
detail bv lKeiiermann et all dl998ll2004h and lLister & Homanl 
d2005l) . for each epoch, the data were calibrated using AIPS 
and imaged with Difmap. The location of each bright feature 
was fitted in the visibility plane. 

We have supplemented these multiple epoch images using 
VLBA 2 cm archive data from observations made at three 
epochs in 2000 (22 January, 8 May, 30 December). These 
later observations included one VLA antenna and were made 
with full tracks in hour angle each lasting about 10 hours us- 
ing 2-bit recording at a 256 Mbps data rate, allowing us to 
reach the thermal noise level of less than 70 /iJy beam" 1 . 

As a check on the robustness of the weak features seen in 
these high dynamic range images, different authors indepen- 
dently made images using Difmap and AIPS respectively and 
obtained very similar results. In addition, to verify the real- 
ity of faint "counter-feature" reported in $3] one of us gen- 
erated three artificial datasets of M87. These model datasets 
were based on the structure of M87 which was observed on 
2000 May 8. They used the identical (M,v)-coverage of that 
epoch and contained a comparable level of random noise. One 
of the models contained an extended counter-feature totaling 
nearly 40 mJy, and two did not. All three model datasets 
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FIG . 1 . — VLBA 2 cm logarithmic false-color image of the M87 jet. The beam, FWHM of 0.6 by 1.3 mas at P.A. =—11° is shown in the lower left hand corner. 
The peak intensity is 1.00 JybeairT 1 and the off-source rms noise 64 /^Jy beam - '. The corresponding dynamic range is better than 15,000 to 1. This image was 
obtained from a full track in hour angle using the VLBA along with one VLA antenna on 2000 May 8. The faint counter-jet is seen extending to the southeast. 



were independently self-calibrated and imaged by another one 
of us without any prior communication about which models 
contained the counter-feature. The image maker was able to 
correctly identify the data with the counter-feature and cor- 
rectly determined its flux density and extended structure. No 
counter-feature was found in the derived image down to the 
noise level from the other two datasets without the counter- 
feature. Given the results of these trials, and the presence 
of the counter-feature in all three observed full-track VLBA 
images as well as in the 20 lower sensetivity 2 cm Survey / 
MOJAVE images , we are confident that the observed counter- 
feature is real and that there are no significant artifacts in our 
image created by the imaging process. 

3. JET STRUCTURE 

In Figure [T] we show the 2 cm image constructed from ob- 
servations made with the VLBA plus 1 VLA antenna. A ta- 
pered image made with twice the beam size from the same 
data shows structure out to nearly 0.2 arcsec. As seen in 
Figure Q] the jet appears bifurcated, starting at about 5 mas 
(0.4 pc) from the core, characteristic of a single limb bright- 
ened cylindrical or conical jet. The M87 jet appears to be 
highly collimated, with re-collimation observed between 2 pc 
where the opening angle is about 16°, and 12 pc where the 
opening angle is only 6° to 7°. As discussed in §|2] FigureQ] 
also shows the existence of weak structure extending away 
from the bright core toward the southe ast. Indications of this 
counter feature were first suggested bv lLv et ail (|2004) on the 
basis of their 7 mm VLBA observations. This counter-feature 
also appears clearly bifurcated but is weaker than the main jet 
by a factor of 10 to 15 close to the core and is only traced out 
to 3.1 mas from the core. Between 3.1 and 6 mas away, this 
feature is weaker than the jet by at least a factor of 200. 

For several reasons we believe that the eastern extension 
may be the counter-jet. Based on their higher resolution 7 mm 



image, Lv et al. (2007) have also detected the counter feature 
and have argued that the true base of the jet cannot be off- 
set by more than 2 mas from the bright core; whereas we 
have detected the counter-feature to be at least 3.1 mas long. 
Also, we note strong circular polarization at a fractional level 
of —0. 5% ±0.1% was detected at 2 cm by iHoman & Listerl 
(2006) coincident with the flux density peak, suggesting that 
this region of the jet has an optical depth near unity (e.g., 
iJonesI 119881) which is characteristic of jet cores. Keeping 
in mind th at there is possib l e amb iguity in the VLBI image 
alignment, Zavala & Tavlor (2003) have measured the spec- 
tral index between 8 GHz and 12 GHz for the region with 
the peak intensity to be a = 0.24 (S ~ v a ) consistent with an 
optically thick nucleus. Since the surface brightness of the 
eastern feature is more than 200 times fainter than the one of 
the brightest feature which we have identified with the core, it 
seems unlikely that it could be the actual core. However, con- 
sidering that no other radio jet has been observed with com- 
parable linear resolution, sensitivity, and dynamic range, we 
cannot exclude the possibility that we are seeing the detailed 
structure of the optically thick core of the jet, and that the 
counter-jet itself is not detected. 

4. JET KINEMATICS 

HST obse rvations of the M8 7 jet made between 1994 
and 1998 by iBiretta et al.l d 19991) have suggested superlumi- 
nal speeds up to about 6c in a region approximately 6 arc- 
sec (0.5 kpc) downstream of the central AGN, consistent 
with what might be expected from Doppler boosting of this 
highly asymmetric jet. The feature, known as HST-1, is lo- 
cated 0.86 arcsec from the base of the jet, and was more re- 
cently observed by Chandra and HS T to vary dramatically at 
both X-ray an d optical wavelengths (lHarris et al.ll2003L 120061 : 
ICheung et alj|2007l) . It was found to have apparent superlu- 
minal speeds up to 4.3 ±0.7c in 20 cm VLBA experiment 
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TABLE 1 
M87 JET COMPONENTS SPEED 



Component 


<- A > 




1/ .> 


Mr 


a 

■ 'app 




(mas) 


(mas) 


(deg) 


Oasyr -1 ) 




(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


C7 


0.43 


0.1 


111.2 


40±6 


0.010 ±0.001 


C6 


0.43 


0.1 


279.2 


13 ±6 


0.003 ±0.001 


C5 


1.4 


0.2 


283.8 


59±11 


0.015 ±0.003 


C4 


2.8 


0.3 


283.5 


104 ±17 


0.026 ± 0.004 


C3 


6.4 


0.5 


282.8 


82 ±26 


0.021 ±0.007 


C2 


13 


1.3 


286.6 


216=1=81 


0.05 ±0.02 


CI 


21 


1.1 


290.0 


12 ±76 


0.00 ±0.02 



NOTE. — Col. (1): Component identifier, Col. (2): Mean radial position relative 
to the core, Col. (3): Typical radial position error calculated from the scatter of points 
from the linear fit (Figure[2), Col. (4): Mean position angle relative to the core, Col. (5): 
Angular radial speed and la uncertainty. Positive speed means that a component is 
receding from the core, negative — approaching the core, Col. (6): Radial speed in units 
ofthespeedoflightand lcr uncertainty. 



by Cheung et al] d2007l) . However, wi t hin 1 arc sec of the nu- 
cleus (80 pc), using HST, Bire tta et~aT] (fl999) measured much 
slower speeds between 0. 6c and 0.8c. 

At radio wavelengths, Bir etta et al] i 19951) reported VLA 
observations which showed a typical speed of about 0.5c 
but ranging up to 2.5c at distances up to about 20 arcsec 
(1.6 kpc) f rom the nucleus. However, previous VLBA ob- 
servations (lBiretta& Junorlll995t lJunor & Biretta! 1 19951) and 
VSOP (VLBA to HALCA) observations dDodson et al.ll2006l) 
found no evidence of mo tions within 5 pc of the core, al- 
though iReid et al] (Q~989) reported an observed velocity of 
0.28c ± 0.08c for a feature located about 20 mas (1.6 pc) 
from the n ucleus. More re cently, based on just two of five to- 
tal epochs, iLv et al] d2007l) reported possible speeds between 
0.25c and 0.4c for three transient features located about 3 mas 
(0.24 pc) from the nucleus. 

We find no evidence for motions faster than 0.07c within 
the inner 20 mas (1.6 pc). In Figure|2]we plot the location of 
seven separate jet features as a function of time, and we sum- 
marize their speeds in Table Q] All the jet components were 
fit with circular Gaussians. The fastest jet speed observed is 
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FIG. 2. — The location of seven separate features which were seen at mul- 
tiple epochs are plotted as a function of time. Component C7 refers to the 
counter-jet. Typical errors of the components' separation are presented in 
Table[T] The large scatter in the positions of the outer two components, CI 
and C2, is due primarily to their large angular extent which makes location of 
their centroids uncertain in our snap-shot images of the limb-brightened jet. 



only 0.026 ± 0.004c, while the counter-jet feature apparently 
moves outward at 0.010 ± 0.001c. The other five features ap- 
pear essentially stationary over the twelve years that they have 
been observed, with nominal upper limits to their speed of 
about 0.07c. Each feature was detected at 17 to 24 separate 
epochs. 

Given that much faster speeds, up to 6c, have been observed 
in the arcsecond scale jet, it is important to consider the pos- 
sibility of temporal under-sampling in our kinematic analysis 
of the milli-arcsecond scale jet, where 6c would correspond 
to ~ 24 mas yr" 1 . Our VLBA sampling interval ranged from 
2 to 8 months, with a median of 5 months. A newly ejected, 
fast component could therefore conceivably travel anywhere 
between 2 to 16 mas downstream in the time between suc- 
cessive epochs, possibly leading to an incorrect component 
cross-identification in Figure [2] However, the large num- 
ber of jet components present at each epoch implies that if 
they have very fast speeds, new components would have to 
be emitted very frequently from the core region, at a rate of 
~ 2 per year. But, the flux density and brightness temperature 
of the core show no indication of such activity, but instead 
vary relatively smoothly over time. Second, there are clear 
gaps in jet brightness between 3 and 5 mas and between 6 and 
10 mas (see Figure|2]and the movie in the electronic edition of 
the Journal) in which jet components are never present at any 
epoch — a highly unlikely situation given a high component 
ejection rate. Third, the brightness temperatures of the indi- 
vidual components CI to C7 are consistent across the epochs, 
which would not be expected if our cross-identifications were 
incorrect. Finally, our 2 cm VLBA movie of the jet 6 , which 
relies on simple linear interpolation across the epochs, shows 
no sudden jumps that would be indicative of temporal under- 
sampling. 

We believe that this interpretation is more robust than the 
results of other radio measurements which typically used only 
two or three epochs. We cannot rule out the possibility of a 
smooth relativistic flow with no propagating shocks or other 
patterns. 

5. DISCUSSION 

Considering that the Doppler boosted relativistic jet of M87 
appears one-sided out to kiloparsec scales, the absence of any 
clear motions in the inner jet is somewhat surprising. We 
measure a jet to counter-jet flux density ratio of between 10 
and 15 (in the year 2000) at a position between 0.5 and 3.1 
mas from the core, and a lower limit of 200 between 3.1 to 
6.0 mas away. Assuming that the jet is intrinsically bidirec- 
tional and symmetric, and that the bulk velocity flow is the 
same as the pattern motion, our maximum reliable observed 
speed of /3 app = 0.026 ±0.004 would imply a jet to counter- 
jet ratio near unit y for the commonly accepted view ing an- 
gle of - 40° (e.g. bwen et 10(19891: IReid et al.lfl989l) . If the 
jet and counter-jet are intrinsically symmetric and oriented at 
an angle of 30° to 40° to the line of sight, the observed jet 
to counter-jet flux density ratios imply that the intrinsic flow 
speed, (3, is 0.5 to 0.6 between and 3.1 mas from the core 
and increases to >0.9 beyond 3.1 mas. Also, we note that the 
rather small ratio of observed jet to counter-jet speeds (Ta- 

6 A movie of the jet covering the period 1999 to 2007 can be seen 
in electronic edition of the Astrophysical Journal. The maps prior to 
1999 have a higher noise, they were not included in the movie. An 
up to date version of the movie, as well as a journal of the observa- 
tions, i mage and visibility FITS file s can be found in our web database 
http://www.physics.purdue.edu/MOJAVE/sourcepages/1228+126.shtml 
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ble[T|i of less than 2.5 is not consistent with simple relativis- 
tic beaming models. We conclude that either the inner jet is 
intrinsically asymmetric, or there are no detectable moving 
features within a rapidly flowing plasma. 

Apparent limb brightening is predicted from analy t ic and 
numerical modeling of relativistic jets (lAloy et alJ 120001 : 
iPerucho et alj [2007) and can be reproduced via Kelvin- 
Helmholtz instabilit y, which is seen in e xtragalactic jets both 
at kilopa rsec (M87, Lobanov et al. 2003) and parsec scales 
(3C273, lLobanov & Zens us 2001). Limb brightening can 
be particularly pronounced whe n the jet opening angle is 
greater than the beaming angle (Gopa l-Krishna et all 120061) 
and especially if the r e is a velocity gradient across the jet 
(Gopal-Krishna et alj 120071) . A two layer "spine-sheath" 
model has been suggested to explain the existence of observed 
strong TeV emissio n from BL Lacs wit h apparently slow 
movi n g radio jets fC hiaberg e et alj l2000t iP iner & Edwards 
I2004t iGiroletti et all 120041) IStawarz & Ostrowskil d2002l) 
and [Ghisellini et al. (2005]) have considered a two compo- 
nent jet having a fast spine which produces the gamma 
ray emission by inverse Compton scattering of the radio 
photons from a surrounding slow layer or sheath. Spine - 
sheath mod e ls hav e also been disc u ssed b y iBridld (I1996T). 
Swain etal.1 d!998l). lAttridge et alj d!999l) . lLaing & Bridlel 
d2002l) . and lCohen etal.N 2007). The central gap seen in our 
VLBA images of M87 , as well as observed TeV emission 
dAharonian et alJ [20061 and as discussed above) combined 
with the lack of measurable motion within the inner 1.6 pc 
(® appear to support this two component model. 

The recently repo rted detection of stron g variable TeV 
emission from M87 dAharonian et al. 2006) also presents a 
problem, since unless the radiating plasma has a large Doppler 
factor, energy losses due to 7-7 pair production will ex- 
tingu ish the gamma ray emission (e.g.. iDondi & Ghisellinil 
1995). Although we do not find any evidence for a fast 
moving jet in M87 close to the central engine, the observed 
TeV emission can be explained in terms of a dual layer 
model with a fast inn er jet and a slower moving outer layer 
(Ghis eTiTni et al.ll2~005l) . In this picture, the inner jet is beamed 
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away from us and is thus not seen in our VLBA images, and 
we only observe the slower outer layer. However, even the 
slow outer layer must move at at least (3 > 0.5-0.8 to be con- 
sistent with the jet to counter-jet flux density ratio as discussed 
e arlier in this se c tion. 

ICheung et alJ (l2007h argue that the superluminal feature 
HST-1, located at a de-projected distance of more than 100 pc 
down the jet, may be the source of TeV emission in M87, 
rather than the base of the jet. By analogy, they suggest that 
blazar activity more broadly may not, as is widely assumed, 
be located close to the central engine. If this is the case, 
we might expect to observe small-scale structure in HST-1 at 
2 cm, but we do not find any milliarcsecond structure stronger 
than 0.75 mJy in that region in the 2000 full-track VLBA data; 
however our data were taken five years be fore the HST-1 flar- 
ing event in 2005. lAharonian et alJ d2006l) commented that an 
origin of the TeV emission in HST-1 would imply an unreal- 
istically small opening angle for the energy source, assuming 
the source was located at the base of the jet. 
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